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Abstract 

We estimate locations of the regions of the percolation and of the 
non-percolation in the plane (A,/?): the Poisson rate - the inverse 
temperature, for interacted particle systems in finite dimension Eu- 
clidean spaces. Our results about the percolation and about the non- 
percolation are obtained under different assumptions. The intersection 
of two groups of the assumptions reduces the results to two dimension 
Euclidean space, R^, and to a potential function of the interactions 
having a hard core. 

The technics for the percolation proof is based on a contour method 
which is applied to a discretization of the Euclidean space. The tech- 
nics for the non-percolation proof is based on the coupling of the Gibbs 
field with a branching process. 

Keywords: Non-ideal gas, Poisson point process, Boolean percolation. 
AMS 2000 Subject Classifications: Primary 82B43, 82B26, Secondary 
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1 Introduction 

A rigorous proof of phase transitions for continuous models of the statistical 
mechanics is still an open problem if the interactions between particles are 
described by conventional in physics potential functions. The first and yet to 
this moment the only example of the rigorous proof of the phase transition in 
a continuous model is the result by J.L. Lebowitz, A. Mazel and E. Presutti 
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in [3]. The potential functions in [3] are a pre- limiting version of the mean- 
field interaction, see [2], having a large but finite radius of the interactions, 
and a four-body stabilizing potential function. 

In the present work, we investigate the percolation properties of an in- 
teracting particle ensemble. We describe a phase diagram of the continuous 
system in the plane (A,/3), the Poisson rate - the inverse temperature. The 
interaction is defined by pair potential functions. We do not prove the phase 
transition driven by boundary conditions as in [3]. However we think that 
the transition: the percolation - the non-percolation can be considered as a 
phase transition relatively, for example, to the conductivity of the matter or 
the velocity of the sound propagations. 

The book [S] gives a rather complete picture of the state of the continuum 
percolation theory for the ideal gas from the mathematical point of view. 
Much attention in [5] is drawn to the Boolean percolation problem for the 
Poisson point processes in My . Points of a configuration of the process are 
considered as the centers of closed balls of a random radius [Boolean radius) 
such that the radii corresponding to different points are independent of each 
other (and also independent of the process) and identically distributed. The 
existence of an unbounded connected component in the set composed by the 
union of all random balls means the percolation. The unbounded component 
is called an infinite cluster. One of the main results in [5] which is related to 
the present article is about the existence of a critical value Ac of the rate of the 
Poisson point processes. Namely, the value Ac distinguishes the percolation 
and the non-percolation, where the last means that with the probability 1 
only bounded connected components exist in the union of the balls. It is 
asserted in [5] that with the probability 1 there are no infinite clusters when 
A < Ac and there exists an infinite cluster when A > Ac. 

We consider the same problem but for a non-ideal gas, which is deter- 
mined by some interaction potential function and the Poisson free measure 
with the rate A. A non-percolation condition was studied in the work [10] for 
positive finite range potential functions. The Boolean radius is equal to the 
range of the interaction. A new proof of this result can be found in [TT] . 

Next we give a brief description of our results not concerning the con- 
ditions. We study the case when the potential function takes as positive 
as negative values. The potential function determines a Gibbs measure of 
which the percolation properties we investigate. In this case a new parame- 
ter enters into the game it is temperature T. By a tradition we more often 
use the inverse temperature /5 = ^. The results we present here outline two 
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regions in the plane {X, (3) of the percolation and the non-percolation with 
probability 1 for a Boolean radius i. We do not seek the solution as precise 
as possible. Our aim is to outline the regions such that they have typical 
forms. Namely (see Figured]): 

The region of the non-percolation can be described as follows. 

There exists a density value such that for any X < there exists an 
inverse temperature /3^(A) such that for any (3 < f3^{X) all clusters are finite 
with Gibbs probability 1. 

The region of the percolation can be described as follows. 

For any density X there exists an inverse temperature /3/(A) such that for 
all (3 > l3^{X) there exists an infinite cluster with Gibbs probability 1. There 
exists a density value A/ such that /5/(A) = if X > X'^ 

Our results provide estimates of the parameter regions separating the 
areas of the existence ^4+ and of the non-existence A- of an infinite cluster. 
There exists a region between and A_ where our result does not give the 
answer on the percolation. 

The result shows in particular that for any small density A there exists an 
infinite cluster if temperature is low enough. Another feature is also that the 
non-existence of an infinite cluster may only be at a small density, A < A^. 
This fact is in accord with the result ([5j) of the Boolean non-percolation for 
the ideal gas. 

The conditions under which we prove the percolation and the non-percolation 
are different. We prove the percolation result in only. It is necessary as 
well, that the potential function has an attractive part. However, we do not 
assume a hard core. Our proof of the non-percolation requires the hard-core 
condition. The attractive part of the potential as well as the dimension of 
the Euclidean space are not restrictions for our proof of the non-percolation 
. Besides, all results are proved for a non-random Boolean radius i. 

The technics for the proofs of the existence and of the non-existence of an 
infinite cluster drastically differ. For the existence of an infinite cluster we use 
technics close to the contour methods (see [T]). The non-existence is proved 
by a coupling of the Gibbs state and a branching process. The extinction 
of the branching process leads to the non-existence of infinite clusters. We 
use a branching process with interactions between offsprings in the same 
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generations and between the generations. The hard core condition prevents 
the accumulations of a large offspring amount which can appear because of 
attractive interactions between the offsprings. 

The section |2] contains the definitions, all assumptions and the formula- 
tions of the main results. All proofs are in the section [31 

2 Model and Results 

The configuration space, the potential function and Hamiltonian 

The non-ideal gas model is a pair {Q,ip). Here Q = {uj} is the set of all 
countable subsets in W such that for any bounded V G'R'^ 

#(cunl/)<oo, (2.1) 

where is the number of points in W. uj is the set of points from R'", 

where particles x G o; sit. The set VL is called the set of configurations. We 
use the standard notations for the restrictions on subsets. If is a Borel set 
in W and uo then ujy = uj r\V and Vty is the set of all configurations in 
V . If V n V = and u & Q then uvuv ='■ V cuy. 
The (T-algebra 21 in is generated by the cylinder sets 

Av,n ■■= W ■ #M =n}<zn, (2.2) 

where ^ is a bounded Borel set in . 

The potential function describes the interaction of the particles. We 
consider pair interactions only and assume that ip{x,y) is continuous and 
satisfies the following properties. 

• Translation invariance: for any (x, y) E M." x M*^ and any z G M'^ it 
holds that ip{x + z,y + z) = ip{x, y). 

Therefore we can introduce the function (p{x), xE M'^, by the equality 
ip{x — y) = ip{x,y), which further we denote with the same symbol 

if{x). 

• Isotropy: if B is an orthogonal operator in W then (p{Bx) = (f{x). 
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• There are two reals / > 0, > such that f < d and 

{= oo, if |x| < /, 
>0, ii\x\e[f,d], (2.3) 
< 0, if \x\ e [d, oo). 

Besides there exists a positive monotone decreasing function ijj and 
g > d such that 



ip{x) > —ip{x) for X > g, (2.4) 

and 

/■oo 

/ = / r''"^ip{r)dr < oo. (2.5) 

The condition (12. 5p was proposed in [T2] . 

• Lower boundedness: there exists M > and Xq such that min^. (p{x) = 
ifixo) = -M. 

Hamiltonian is 

Hiuj) =Y,v{x- y) (2.6) 

which describes energy of configuration uj. The above expression is formal 
since the sum does not exists. The energy of uy G with boundary 
condition r G Qy^ is 

H{uv\t) =H{iJv)+F{uv,T) := J2 Vix-y)+ Yl v{x-y).{2.7) 

The last sum in the above expression might be infinite. However, if / > 
and (12.51) holds then J^xeujy yer'^i^ ~ 2/) < oo for any finite V, and any uy 
and T. 
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The reference and the Gibbs measures 



The reference measure 11 is defined as Poisson one on [Q, 21) with intensity 
A > : 

n(^v,n) = ^^e-^l^l, (2.8) 

where \V\ is the volume of V (see f l2.2p ). The Gibbs measure P^'^ on (17,21) 
is determined by the Gibbs reconstruction method of the reference measure 
(see m). 

To define P^'-** we introduce a Gibbs specification 

which is a family of the Gibbs reconstruction of the measure 11 in finite 
volumes V given a conditional configuration r and the inverse temperature 
(3 G M+. The measure Py^ has the following density Py'^ with respect to the 
measure 11 : 

ftA/ X exp{-l3H{ujy I r)} 

Jnv ew{-pH{uj I r)}n(d^y) 

It is assumed that the densities py^i^uy) are defined for boundary configura- 
tions r such that (12. 7p is finite. 

We assume some conditions for the existence of the integral in (12.91) and 
for the existence of at least one of Gibbs measures P^'^ corresponding to the 
specification (12.91) (see [7] or [B]). If H{-\t) is finite not for all r then the 
existence conditions is such that the Gibbs measure P^'-*' is concentrated on 
a set of the configurations, where H{-\t) is finite, when r from this set (see 

Further we use the notation Hy{-\-) for the energy of configurations from 
Qy with a boundary condition. We shall often omit some indices and shall 
write Py instead of Py^ and P instead of P'^'^. 



The percolation 

Any ordered sequence of particles from a gas configuration uj we shall call a 
path. We say that two points x,y eM'^ ^-percolate with respect to u if there 
exists some finite path n = {xi, X2, ■ ■ ■ ,Xn} Co; such that \xi — < i for 
alH = 2, . . . , n, and \x — xi\ < i, |y — < £. The set vr is called a i- cluster 
in UJ or simply a cluster when i and u are fixed. We say that x,y in a 
cluster TT if they ^-percolate with respect to u. We shall denote by (x y) 
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the event that x and y are in a cluster. If there exists an infinite cluster 
starting at x G M"^, we denote this event by (x oo). The probability of 
the event (0 oo) we call i -percolation function or percolation function and 
denote 9i{P, A) : 

ei{p,\) := P'^'\0^^ oo). (2.10) 
The main result on the non-percolation 

The result on the non-percolation is proved for / > (the hard-core condi- 
tion) and any u (any dimension of the Euclidean space). 

Theorem 2.1. For any i > f there exists \J , < \J < oo, and a function 
Pi{X) defined on the interval (0, A^) such that 

1. < < oo on the interval (0, A^), 

2. /3^{X) t oo as A i 0. 
Let 

A_ = {iX,/3): A<A7, /?< /37(A)}. 

Then all clusters are finite with the probability 1 if {X, (3) G A_, that is 
6i{l3,X) = for those (A,/5) (see FigU^. Moreover the expectation of the 
cluster size is finite. 

The main result on the percolation 

The next theorem is proved for the case u = 2. Now the Boolean radius is 
bounded below, however the hard core is not necessary. 

Theorem 2.2. For any i > 2^/2d there exist A^, < A^ < oo, and a 
function /9/(A) defined on the interval (0, A^) such that 

1. < j3^{X) < oo on the interval (0, A/), 

2. Pl{X) t oo as A i 0. 
Let 

A+ = {{X,(3): A<A+,/3>/3+(A)}U{(A,/3): A>A+,/3>0}. 

Then with the probability 1 there exists an infinite cluster if {X, P) G that 
is 9i{p, A) > for those (A,/3) (see FigUD- 
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Figure 1: Percolation and non-percolation regions. 



The next corollary joins the results of theorems 12.11 and 12.21 However, 
the claims of the corollary hold under the most restrictive assumptions from 
the assumptions of both theorems. 

Corollary 2.3. Let u = 2 and / > 0. For all i > 2\/2d there exist two 
positive reals Xj and A/, and two functions: defined on (0, A^), and 

/5/(A) defined on (0,A/), such that 

1. < A7 < A^ < oo, 

^. < (3'^{\) < 00 and < /3/(A) < 00 on their intervals of the defini- 



tions 



3. 



/?7(A) </3+(A) on (0,A7), 
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Pi (A) t 00 and /5/(A) t 00 as A | 0. 



Let 



A 



{(A,/5): A< A7, /3< /57(A)} 



and 



A+ = {{X,P): A<A+ /3>/?+(A)}U{(A,/3): A > A+, /5 > 0}. 



Then 
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1. with the probability 1 all clusters are finite if (A,/3) G A_, 

2. with the probability 1 there exists an infinite cluster if (A,/5) G Aj^, 
(see Fig\Ti). 

Moreover the expectation of the cluster size is finite if (A,/3) G A_. 



3.1 The proof of the non-percolation 

It is essential for the proof that the potential function has the hard core. 
However, the arguments in this subsection do not depend on the dimension 
of the space. 

The next lemma is an obvious consequence of the hard core condition and 
the inequality fl2.5p . 



Lemma 3.1. Let V (^My be a finite volume and let a be a configuration in 
V. Then 



We obtain the non-percolation result by a coupling of the Gibbs measure 
P^'^ with a branching cluster process. The method was described in the work 
[H]. Here we use the analogous idea. 

Note that due to the hard core condition, one considers the ^-percolation 
only when i > f. Informally the idea of the proof is the following. Suppose 
that i? C u; is a ^-cluster with #(-R) > 1. Let us choose some particle xo 
from this cluster, xq G R. Let i?^^^ C i? be the set of all particles from R 
such that the distance between xq and any point from R^^'' is less than or 
equal to i. The set R^^^ is not empty because #(-R) > 1. Next we choose a 
set R^^\ where R^"^^ C R \ (R^^^ U {xq}), and which is the set of the particles 
such that for any v from R^'^^ there exists at least one point w from R^^^ at 
the distance no greater than i. We can call v an offspring of w. If the set 
R \ U U {xo}) is not empty we can choose a subset R^^^ with the 
similar properties, ect. Iterating the procedure we will obtain the following 
representation of the cluster R : R = U°ZqR^^^ (here R^^^ = {xq}). The set 
R^^^ we call i-th generation. The set R^"^'^^ generates the set R^"'\ 

This branching construction brings us the idea of a branching process, but 
there are two peculiarities that differ our process from the ordinary branching 



3 Proofs 




(3.1) 
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process. First, note that it is possible for one offspring to have different 
parents. Thus we do not have here a branching tree and it means that the 
independence of the offsprings does not hold. Second, keeping in mind the 
coupling, we define the transition probabilities of the branching process by 
Gibbs measure Py- Thus for £ sufficiently small it is possible that some 
generation R^^^ interacts with precedent generations R'^^\ k < n. 

Next we describe a rigorous construction of the cluster branching process. 

3.1.1 The cluster branching process 

We describe a path of the process and its distribution P. 

Let xq be some point from W^. We shall construct a sequence which 
describes the generations. Together with the sequence (i?*-"^) we define the 
sequence where = Ui=oR^'^ and the sequence of the occupied areas 

Bn = U^(zE„_iB£{v). The set we call an environment. 

Initial step. Let R^^^ = {xq} and Eq = {xq}, Bq = 0. 

First step. Let i?^^'' = {x^i \ . . . ,x^k^} be a set of the particles in the ball 
Bi{xq) with the center xq and with the radius ^. The set R^^^ is the offspring 
set of Xq. Then Ei = EqU R^^^ and Bi = Bi{xq). Bi is the occupied area by 
the offsprings of xq. No particles of further embranchments appear in Bi. We 
define a conditional probability density p of the measure P with respect to 
the same Poisson measure 11. The density with respect to 11 of the offspring 
set R^^^ of the ancestor xq is 

p{x^^\ . . . ,4;) I Eo) = ^ / e-^^(^^-)n(do;), (3.2) 

1 

where Qb^ is the set of all configurations where particles "live" outside of the 
ball Bi, and 

Z{Eo) = [ e-^^(^°^")n(da;). 
Jq 

We use Gibbs measure P^'^ for the definition of p. In fact, all defined prob- 
abilities and further calculations assume a big volume V, where we consider 
all configurations. It means that we use the measure Py instead pl^'-'^. For 
example, in (13.21) the integration is taken over Qb^hv- Therefore p depends 
on V. However, in what follows all estimates do not depend on V, and hence 
can be considered as the estimates in the infinite volume. We do not mention 
the volume V in the further calculations except cases when it is required. 
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Using p we can calculate, for example, the probability to have k offsprings 
of Xq : 

P(#(i?«) = A; I Eo) = / I Eo)n(di?(i 

_ ^ l^ll ^-AlBil / ^/ (1) ^(1) 



k\ 



Second step. In this step we describe the embranchments of all particles 
from R^^\ Any particle branches according to some order introduced in 
R^^\ We shall construct the set i?^^^ in according to the chosen order. Let 
^(2,1) ^ . . . ,x^^^^^} C Bi{x^^^) \ Bi be the set of offsprings of 

The offsprings of x^i^ cannot be situated in Bi since Bi is occupied by the 
offsprings of xq. Let E^2,i) = Ei U i?^^.!) and S(2,i) = Bi U Be{x^^^). The 
probability density of the offsprings of x^^^ is 

1 



p{xf^'\ xgf I E,) = ^ / e-^^(«'"^^^^^^-)n(d^), (3.3) 



'(2,1)^ 



where 



Z{Ei) = I e-^^(^^^")n(du;) (3.4) 

^1 

(see (O)). 

Suppose now that k < ki particles from R^^^ are already branched and 
R^'^'^\ . . . , is a sequence of their offsprings. Hence we have the environ- 
ment (all already living particles) -E(2,fc) = EiU |Ji=i -R*-^'*^ and the occupied 
area -B(2,a:) = -Bi U IJi=i -^^(^!^^)- Note that i?(2,A:) is the ^-neighborhood of 
EoU{x^^\...,x^^}. 

Let now the next point x^^]_-^ be branching. Let 

be the set of offsprings of x\!'^-^. Then E{2,k+i) = E{2,k)^R^'^'^^^^ and -B(2,fc+i) = 
-8(2, fc) U -B£(x[,^_^;^). The probability density of the offsprings is 

p (xf^'-'\ . . . , xg:::^ i e,.,)) = / e-/^-(-(--)v-)n(do.). 

(3.5) 



11 



We obtain the next generation R^'^^ = Ijfi]^ i?'-^'*^ after the embranchments 
of all points from R^^\ Let E2 = -E'(2,fci) and B2 = B^2,ki)- 

{n + l)-th step. To construct R^"'~^^^ from R^^^ we follow the same scheme. 
Let i?*^"-* = {x^^\ . . . ,x^i^]}- The particles from i?^"-* are branching according 
to some chosen order in R^"'\ Suppose that k {k < k^) particles from R^"^^ 
are branched. Hence we have the sets (where i < k), E(^n+i,k) and 

B{n+l,k)- 

Now, let 

be the set of offsprings of the branching particle a;^*^]^. Then E(^n+i,k+i) = 

E(„+i,fc) U i^C^+i-'^+i) and B^n+i,k+i) = B^n+i,k) U 5^(4'?i)- The probabilistic 
density is 

/ (n+l,fc+l) (n+l,fc+l) I p 

g-/3H(i?(„+i,,+i)Va;)n^^^^)_ (3.6) 



Z{E^n+l,k)) JufB, 



The above iterative steps describe a path and the transition probabilities 
of the cluster growth process. Note that the offsprings can depend not only 
on the preceding generation, but on all previous generations. 

The question we are interested is when the cluster branching process 
extincts with the probability 1. It is known that the extinction condition 
of the ordinary branching processes is formulated in the term of the mean 
number of the offsprings. We can expect that the extinction of the cluster 
branching process is controlled by the offspring mean value, as well. We 
show further that if the mean number of the offsprings is uniformly less than 
1 over all possible cluster configurations of the previous generations, then the 
cluster branching process will extinct with the probability 1. 

In the next lemma we prove that there exists a region in the plane (A, (3) 
such that the mean offspring number of one ancestor is less than 1. Then, 
we show that this condition is sufficient for the extinction of the cluster 
branching process. Moreover we show that the mean value of paths of the 
cluster branching process is finite. 
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Lemma 3.2. There exists and a function (3 (A) such that for any A < 
and (3 < P~{X) the expected number of the offsprings ^(i?'^"'^'') of the ancestor 
less than 1, uniformly overn,k and over the environment E(^n,k-i)- 

sup < 1. (3.7) 

n,k 

Proof We give an estimate of the probability the point x^" to have 
exactly K offsprings, that is = K. By the definition of the offspring 

density (13. 6p we have to estimate the following integral. Let : = 

Bt{x^k~^^) \ 5(„,fc-i) then 

P(#(i?^"''^)=i^ I ^(n„fc-l)) (3.8) 

= Ic , p(i?^"'') I i?(„,fc-i))n(di?("''=)) 

We shorten some notations in the further calculations. Let R := and 
E := ^\ Than the integral (13.81) can be represented as 

P(#(i?("='=)) = K I i?(„,fc_i)) = P(#(^) =K\E) (3.9) 

= [ [ e-^^(^^^^-)n(da;)n(d^), 

where 

Z{E) = I e-^^(^^")n(du;). (3.10) 

Hamiltonian in (13. 9p can be represented as 

H{R yEVLu) = H{R) + H{E V cj) + F{% EVlu). (3.11) 
Let Cm be a sequence of the strips 

Cm = {x G M'^ : m < |x — x^" """^1 < m + 1}, 
where m G Z_|_ and m > g, and 

Co = {x G M'^ : £ < |x - x^:""^^! < mo}, 
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where rriQ = min{m : m > g}. The volume of Cm is 

\Cm\ = KP^_i{m), 

where P,y_i(r) is a polynomial having its power — 1, and the leading coef- 
ficient is u — 1. Let Qm = ^Cm\B{„fc)- Let uj G flm then the number of the 
particles of the configuration u \/ E in Cm has the following bound 



# (c^VE)c„J < 



IC 



ay 



< 



where {uj\/ E)c^ is the restriction of E to Cm- We have then the following 
lower bound for the interaction energy of {u V E)c^ and R consisting K 
offsprings, = K, 



F({ujyE)c„^,R)>-i,{m)K 



.2^P._i(m) 



The energy of the interaction oi ujM E and R is estimated as 

F (ujy E,R] > - KMno - -— ^ P^_i(m)V^(m) 



> 



> 



m=mo 

KMno —I Pu-i{r)tp{r)dr 



mo 



K{Mno + —I 



where no 



ICol 



(see (123])). Let ni = Muq + j^L 



By ([31]) and fim]) we have 
P{R : #(P) = K I E) 



< 



-I3H(R) 



Z{E) J [nm , ■.#{r)=k] 



n(dP) 



-l3H{Eyu}) 



n(dcj). 



{n,fe) 



The cluster energy H{R) we can estimate very roughly as 



K 



-H{R) < M( ] < MK^ 



(3.12) 
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Thus 



P(i? : = K\E) 



< 



< 



Z{E) ijiJen . („_i)y#{ij)=x} 



n(di?) 



-f3H{EVu}) 



(n,k) 



^l3{Km+MK^) 



< e 



Z{E) K\ 



-l3H{Eyui) 



n(dcj) 



n(da;) 



(n,k) 



l3{Kni+MK^) 



K\ 



(3.13) 



In the last inequahty we use 

1 



Z{E) Jn 
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(n,k) 



Besides in the last inequality in fl3.13p we estimate the volume \B{x 



(n-l)x 



by 



\Bf(x 



(n-l)^ 



1. Note that B(x 



(n-l)^ 



can be empty, that means that the particle 
x^j^^~^' has no offsprings. It leads to zero of the probability we are estimating. 
We do not use this possibility. 

The right hand side of (13.131) does not depend on the volume V. 
The maximal number of the particles in the ball Bi{x^^ is ub = 
|S^(x[."~^'')|/(k(//2)'^). Thus using estimation (I3.13p we can estimate the 
mean number of the offsprings: 

E(#(^)|E) < J^ke^f^n.+P'^^i ^MrnXk — 



fc=0 



Let 



/5,-(A) 



4lnA-^A-ln 
A A 



~A~ 



(3.14) 
(3.15) 



where 

A = nB{MnB + ni). 

The function is considered on the interval (0, A^], where A^ is a root 

of the equation 

- In A - Kt\ - \n{Kt) = 0. (3.16) 
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When A < and P < P (A) we obtain 

E(#(^) I ^) < 1 (3.17) 

uniformly over the environment E. □ 

Next we show that if the mean number of the offsprings is less than 1 (see 
fl3.17p ) then the mean size of the cluster is finite. Indeed, if f3 < f3~{X), then 
there exists e > depending on P and A such that E(#(i?) | < 1 — e, and 

oo 

E(#(i?("))) = E(5^/^#(^(„-.))=,j#(i?W)) 

k=l 

oo k 

k=l i=l 

oo k 

= ^ Ef J|.^(^(n-l)) = fc|E(#(/g^"''^)|E(„^j_l)) 

k=l i=l 

oo k 

^ (1 - ^) E E E(/{#(/?("-))=.}) = (1 - 6)E(#(i?("-i))). 

k=l i=l 

It means that 

E(#(i?("))) < (1 - e)", 
and we see that the mean cluster size is finite, 

oo oo 

e(5^ #(/?("))) = ^E(#(/2W)) < 1/6. (3.18) 

n=l 71=1 

Since the path i?*^"^ is a ^-connected set of points in V and the estimate (13.181) 
does not depend on V, the probability of infinite clusters is 0. 

This proves that the cluster branching process is degenerated. 

Any path of the cluster branching process starting from xq is 

where is the number of the generations. is finite since (13.181) . The 
relation (I3.18P can be rewritten as 

E(#(i?)) = y2k [ p{E\ EoMdE) < oo. (3.19) 

. J{E: #{E)=k} 
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The coupling 



We explain next the coupling of the Gibbs field having the distribution P^^'-^ 
and the branching cluster process having the distribution P. To this end we 
represent any finite ^-cluster of a configuration u; as a path of the branching 
cluster process. 

Let 7 C cuq G f2 be a finite ^-cluster in a configuration uq. It means that 
Bi/2{l) = UxejBi/2{x) is connected component in Bi/2{u;o) = UxeujoBe/2{x). 

Choose a particle xq G 7. We can consider the configuration 7 as a 
branching process path E = (En), En C En+i, starting at the particle 
Xq {Eq = {xq}), and such that 7 = UEn- The construction was described 
in the section 13.1.11 By the described iteration we obtain also a sequence of 
generations (i?'-"^) such that 7 = U^i?*^"-'. Any generation R^"'^ is a set of the 
offsprings of 

Next we find a relation between the distribution of the cluster branching 
process and the Gibbs measure. Consider the event 



A.f = {Lj: 7 C u;, Bi{j) n 7 = 7}. 



We can define a density x of this event with the respect to 11. 




(3.20) 



The density p(7) of the path 7 is 



TV 



p(7 I Eo) = n P(^^"^ I E, 



'n-l 



(3.21) 



n=l 



It follows from finiteness of 7 that is finite. 
It is not difficult to verify that 




(3.22) 



Then 



x(A,) = p(7|i?o)^ 



(3.23) 



Let (Vm) and (V^) be sequences of the boxes 



Vm = {xeM" : \x\ < m}, Vm = {x eW : |a;| < m - i}. 
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Define the sequence (7^) where 7^ = 7nKn- For any 7^ we have the relation 

m 



Finiteness of the expectation of the ^-clusters follows from fl3.23p . Let 

= {7 C : #(7) = k, xo e 7} 

be the set of all ^-clusters containing the particle xq and having exactly k 
particles. 

Consider the set of the configurations 



and its Gibbs probability 



1 



e 



n(do^)n(d(7V0B,w)), 



where (pB^'j) is the empty configuration in the ^-neighborhood of 7. Hence 

Pv{n^'o) = [ p(7 I xo)n(d(7 V 0B,(7))) 



Irk 

The mean value of the size of the clusters 7 then is 

00 

k=l 

CO „ 

= p(7 I a;o)n(d(7 V < 00 

k=i 

since (13.191) . This proves Theorem 12. II □ 
3.2 The proof of the percolation 

The proof of the existence of infinite clusters is based on technics which is 
close to the contour method in the lattice models. To apply the method 
we discretize splitting it into squares. A c-contour around is a set of 
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empty (without particles) squares surrounding 0. The main fact we prove is 
that the probabihty of a c-contour decreases exponentially with its length. 
It leads to the finiteness of the number of the contours surrounding 0. 

In our proof of the percolation, the essential assumption is that the space 
is two-dimension. The hard core is not used. 

Divide into square cells S = {Sl^j^i^} of the linear size equal to q. 
Suppose that the left-lower corner of any cell 5'^^^. has coordinate {kq,lq), 
where {k, I) ^7? . So we have a natural order of the cells. The point ci^k,i) = 
i~2^Qi ^^?) is called the center of the cell S'^''^ Two cells S"^^^ and S^j^, 
are neighbors if either k = k' ± 1 and / = /' or / = /' ± 1 and k = k'. Let 
(c, c') be the line connecting the centers c = C(^k,i) and c' = C(^k',i') if ;) and 
Slj^i are neighbors. Let V = {S^^j^ ^-j} C 5 be a finite subset of the cells and 
^i'P) = {c(fc,o : S^^j^, i,-^ G V} be the set of all centers of the cells from V. For 
every set V of the cells we consider the graph 

Gv = (c(P),r(P) = {{c,c') : c,c' G C(P)}) 

having C{V) as the vertex set and r('P) as the bond set of all bonds con- 
necting neighboring cells from V. A set of cells V is connected if the graph 
G-p is connected. 

A set of cells TZ is called contour if the bond set r(7^) is homeomorphic to 
the circle. The number n{TZ) of the cells in a contour TZ is called the length 
of the contour. 

If P is a set of cells then W{V) = [jg^-p 5 C M2 js the support of V. 

All contours we consider further surround G M^. Therefore we often 
omit mentioning this. Let a; G be a configuration. If a contour TZ is such 
that ujnW{7l) = then we call it a c-contour with the respect to uj or simply 
a c-contour. 

The proof of Theorem 12.21 is based on the following lemma. Let Q^{71) 
be the event (the set of configurations) such that the contour TZ is the c- 
contour with respect to any u G Q^(7Z), that is fi°(7^) = {ujwcr) = 0)- Let 
Q = Q]v^ny be the set of all configurations out of W(JZ). Then 

^]°(7^) = {uj = (pwiii) V S3 : to en}, 

where (pwin) is the empty configuration in W(7Z). 
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Lemma 3.3. Let the cell size be q = 2d + S, where S is a small positive 
number. There exist constants a G (0, 1) , c(/5) > and G{P, A) such that for 
any h > there exists a function /3//j(A) > defining the domain 

{{X,P): /3>/5+,(A)} 

where the following relations hold: 

1. G{p,\) > h, 

2. for any c- contour TZ 

P^^\n\n)) < c(/3)e-"(^)"^^'^). (3.24) 

The probability that there are no particles in set W(TZ) which is the 
support of the c-contour TZ exponentially decreases with the contour 
length. 

Proof Let y be a volume in containing W{71) = IJ5e7^ ^- order to 
estimate the probability of event f2°(7^) (see Figure [3]- A) we construct a event 
Q^{TZ) by adding particles in the c-contour TZ (see Figure [3]-B). That allows to 
obtain the lower bound of probability Py(f2^(7^)) of the form Pv{^l^(TZ)) > 
e""^(^'^)Py(fi°), where n = niJZ). Substituting the probability Pv{VL^{TZ)) 
by 1, we immediately obtain (13.241) . 

The probability of fi°(7^) 

Recall that VfiJZ) is the event composed of the configurations in V containing 
c-contour TZ. We assume that the boundary configuration out of ^ is r = 0. 
The probability of the event then is 

Pv{QP{TZ)) = ^ [ exp{-pH{uj)}U{du) (3.25) 

-XAn r 

= / exp{-/3i7(^)}n(d^), 

where A = {2d+5Y is the volume of any cell, Vt is the set of all configurations 
iYiV\W{TZ). 

Let (pwin) be the empty configuration in the region WiJZ). Any con- 
figuration oj G VL^iJZ) is the composition of (pwCR.) and a configuration uj in 
V\W(JZ), UJ = 0vi/(7^) VtD. The pre-integral factor in (I3.25P is the integration 
result over (pwin)- 
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Figure 2: The potential function. 
The construction of the event n^{Tl) 

Let m be a positive number such that m < M. For such m there exist positive 
numbers a and e such that e < 6 and 

(p{x) < —m for all |x| G [a, a + e], 

(see Figure [21) 

Let 7 = IJ^^ c')er(7e) ('^' ^^^^ composed of the bonds (c, c') G 

r(7?.). The length of 7 is equal to {2d + 5)n{TZ). We introduce a direction 
on 7. Let it be counterclockwise. The points of any finite subset D of 7 
we numerate along the direction of 7, that is -D = {xq, x^}. The point 
Xj+i G -D is the first point which can be reached from Xj in the direction of 
7- 

Let Ba+^{xi) be the closed disc of the radius a + f with its center at 
Xi E D. If a + I < 4(i + 25 then Ba+£{xi) fl 7 is connected. In the case 
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a + I > 4d + 26 the set Ba+^{xi) fl 7 can be unconnected. Let C{xi) be a 
connected component of fl 7 containing Xi. 

Now we consider D only such that Xj+i G C(xi) and — Xj| = a + | 

for any i < k. Every pair x,y & D such that |x — y| = a+ 1 we call connected. 

The number 77(-D) of the connected pairs in D takes one of two values 



(2(1+5)11(11) 



- 1 < r^{D) < #(/}). 

The number #(-D) of points in any D is not greater than 

Let D = D{n) be such that = sup^ #(5). 

An inverse estimate is in the next 

Lemma 3.4. There exists a > such that for any contour TZ the number 

#(D(7^)) > an{TZ) (3.26) 

^/r^(7^)>^. 

Proof We consider two different cases distinguished by the following 
inequalities 

case 1. a + I < 4(i + 25, 
case 2. a+^>Ad + 25 

The length of the piece of C{xi) C 7 between Xi and Xj+i can be 

estimated as following 



^ {'•^ + 2) ^'^^ ^^^^ 

\2d+& 



s{Xi,Xi+i) < S / , /a+f , A . 

(« + l)2^ + l m the case 2. 



Since 



we can take 



a > 




in the case 1. 
in the case 2. 



□ 
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Let be the disc of the radius | centered at x G -D and U = 

UxeDB^{x). Every disc B^{x) is called a bead and the set U is necklace. 
The set 

Eu = {a enu: Vx e D, 4^{anBs.{x)) = 1, a H t/"^ = 0} 

is a set of configurations all particles of which are located in the beads only, 
one particle in every bead. 

The configuration set Q}{TZ) contains configurations composed by the 
joint of three configurations: 

9}{TZ) = {uji = (jy (t)w{n)\u V 2; : a G Sf/, 25 G } (3.27) 

where (t)w{n)\u is the empty configuration in W{1Z) \ U and u E VL are 
configurations in WiTZy. 

The lower bound for P^-^ 
The probability of Q}{TV) is 

Pv{^\n)) = ^ [ e-^^-(-)e-'^^-(=^)e-'^^(-'^)n(d(aV0vy(7e)wV25)). 

The energy F{a, u) of the interaction of a and uj is negative because the 
distance between any particles of uj and of a is greater than d, hence e~^^^'^'^^ 
is greater than 1. Since (t)w{n)\u = then 

-XnA X^#{D) r r ~ 

Pv{n\n)) > / e-^^^(")n(da) / e-'^^^('^)n(d25). (3.28) 

To estimate Hy{a) remark that there exist at least ^{D) — l the connected 
pairs in D. The interaction energy of any connected pair x,y E D is estimated 
from below as 

(p{x — y) < —m. 

Other pairs x,y E D which are not connected, that is [x — ?/| ^ a+|, interact 
with a negative energy. 
Hence 

J e-'^^^('^)n(d(a)) > e"^'^^*^^)"!) (^J e-^#(^) 
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and it follows from fl3.28p that 

Defining c{P) = e^'" and 

G(/5,A) = /3m + lnA + ln(^^^ (3.29) 

we obtain 

Pvi^'m) < exp{-#(D)G(A A)}exp{/3m} 

= c(/5) exp{-n{n)aG{p, A)}. (3.30) 

The inequality (13.241) holds in the infinite volume since the right hand side 
of ( 13.30p does not depend on V. 
Taking 

/5t.W = -:^MA)-^in (^)+^ 



m m \ IQ J m 

we complete the proof of Lemma 13. 3[ □ 



Next we define 



/3+(A)=/3;^(A), (3.31) 



where c is a combinatorial constant such that the number of the contours of 
the length n surrounding G is not greater than c". It is known that 
c < 3. Let A/ be the solution of the equation /3/(A) = 0. Define the set 

= {{f3, A) : A < A°, /3 > /5+(A)} U {(/5, A) : A > A+, f3>0} 

Lemma 3.5. // {(3, A) G then with the probability 1 there exists only a 
finite number of c- contours surrounding G M^. 

Proof Let Vf'iJZ) be the set of all configurations containing a c-contour IZ, 
the empty contour which surrounds G M^, and let VL^ = U7e- n(7^)=n ^°('^)- 
Then = fio and 

J2P^'\^n) <exp{/5m} J]exp|-n(^aG(/5(A),A)-ln(c))} < oo (3.32) 

n>l n>l 
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if (P, A) G A+. It follows from (^M) that 



P''' n U = 0- (3-33) 

\ m n=m / 

The inequality fl3.33p means that with the probability 1 there exists a finite 
number of the empty contours surrounding 0. □ 
Let u; be a configuration. The set = IJxew-^^/2(^) ^e represented 
as the union of £/2-neighbor hoods of -^-clusters which are connected compo- 
nents. 

We define now a h-contour (Boolean contour). Assume that there exists 
a line L C surrounding G such that Kl (Iuj = 0, where 

KL=[j B,{x). 

The set is called a b- contour surrounding or simply a b- contour. 
The i/2 neighborhood of any ^-cluster does not intersect L. 

Lemma 3.6. Assume that is split into cells of the linear size q. For any 
b- contour K with a radius r,r> V2q, there exists a c- contour TZ such that 
W{n) C K. 

Proof The proof is based on the following simple observation: if we cast 
a coin of the radius r on the plane divided into the square cells S, then 
there exists a cell which will be covered entirely by the coin. Moreover if the 
center of the coin lies on a boundary of two cells or four cells (one point) 
then all those cells are covered by the coin. □ 

We say that two 6-contours are different if c-contours included into them 
are different. Since the number of the c-contours is finite with the probability 
1 then the number of different 6-contours is finite as well. Therefore there 
exists an infinite component in Q^^ for almost all tu. □ 
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7d+d 
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